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ABSTRACT 

Possible applications of the principle of similarity to fluid mechanics is de- 
scribed and illustrated. In correlating thermophysical properties of fluids, the 
similarity principle transcends the traditional corresponding states principle. In 
fluid mechanics the similarity principle is useful in correlating flow processes 
that can be modeled adequately with one independent variable (i.e., one-dimensional 
flows). In this paper we explore the concept of transforming the conservation 
equations by COMBINING similarity principles for thermophysical properties with 
those for fluid flow. We illustrate the usefulness of the procedure by applying 
such a transformation to calculate two phase critical mass flow through a nozzle. 


KEY170RDS 

Fluid mechanics; principles of similarity; thermophysical properties; corresponding 
states; thermodynamic critical point; two phase choked flows. 


INTRODUCTION 


Many dynamic processes can be modeled in terms of one physical coordinate dimension 
and the success of such modeling is well known. However, in doing so one must 
realize that not all boundary conditions associated with higher dimension coordi- 
nate spaces can be satlsifled; Implying the use of a set of suitable constraints, 
a priori, which "average" the dynamic character of the flow field, For example, in 
fluid mechanics the representative velocity, heat flux, and friction factor can be 
written as 


u = f u(Aj^, t) ■ n dA dt; q = j q(Aj^, t) • n dA dt 


^Suppoted in part under NASA Grant NRG- 2 1-002 -344 to the University of 
Maryland , 
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where Aj^ is a surface and a specific time Interval, In order to formulate 

a set of constraints, one postulates the existence of principles of similarity 
between the thermodynamic parameters and the conservation equations, such as con- 
stant or power law dependent fluid properties, or coordinates which behave in such 
a way as to satisfy separation of variables. Often such relations imply a "slug 
flow" velocity profile, or constant heat flux or an overage friction factor based 
on the inlet and exit pressure differences. Under these conditions, only a limited 
number of boundary conditions can possibly be satisfied. 

A family of functions of a number of physical variables is said to satisfy a simi- 
larity principle if the function can be mode identical by reducing the functions 
and variables in terms of suitably chosen scale factors. Thus when such scale fac- 
tors can be found as units to nondimenslonolize (normalize) other variables, the 
profiles plotted in these dimensionless coordinates become congruent, For instance, 
in fluid mechanics when the momentum equation, cited below, is made dimensionless 
using the fluid velocity (u), density (p), viscosity (q) and a characteristic 
length (1), it follows that two flows become similar when their Reynolds numbers 
are equal 


Re^ = 






= Re, 


and this is known as Reynolds principle of similarity. For thermophysical proper 
ties , van der Waals' principle of corresponding states assumes that property sur- 
faces become similar when the variables describing these surfaces ore properly 
scaled, in terms of the critical point parameters. 



where Pj^ = P/P^) % “ T/Tq and Vr = V/Vc* And in ma thematics . similarity solu- 
tions are based on suitable scale factors which stretch the coordinates and reduce 
the governing partial differential equations to ordinary ones. When confronted 
with a complex problem, it is often quite useful to determine the exact solution 
subject to restricted geometric conditions similar to those of the original prob- 
lem. Boundary layer theory is rich in such ,an approach where most similarity solu- 
tions are also one dimensional (Bird and colleagues, I960; Cebeci and Smith, 1974j 
Hansen, 1964; Hinze, 1975; Moore, 1964; Schlichting, 1968). 


The reason for similarity and/or one dimensional modeling can be readily understood 
by considering the conservation equations of fluid mechanics (Bird and colleagues, 
1960; Cebici and Smith, 1974; Hinze, 1975; Schlichting, 1968; Zeleznik, 1976). 


Contunlty: J + _ = 0 


+ pn^u, = + PS^ 

j 


Momentum: ^ + -r— - nu,u, = — 

ot dX. i j dx, 
j i 


3pH g g g 

Energy: — + ^ pu H = ^ (u^T - q.) + ^ 

J J J J 

State: P = P(p,H,x^) 
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where H = 11 + u.u. 

0 2 11 

These equations are difficult to solve without some 8im.plifylng assumptions, even 
for potential flows. With added complexities, such as diabstic flows with real 
fluid properties and severe gradients, the mass and energy fluxes, ttj and qj 
respectively, become most difficult to determine as BOTH T^j and qj are related 
and yet constrained Ivy some boundary. Much of our information for these related 
constraints is constructed from engineering data for the fluids water and air. 

Still our knowledge is limited and we are forced to make simplifying assumptions. 

For the cases where tractable analytical and/or numerical treatments are available, 
regardless of simplifications, we must have accurate, consistent thermophysical 
property information to transform those solutions into physical reality. SometimeG 
this Information is limited or unavailable, yet the success of the entire trans- 
formation depends on its accuracy. So we need to develop methods to predict 
thermophysical properties and techniques to relate these results to an experiment- 
al body of data for fluid a to a similar experiment with fluid p without essen- 
tially reproducing the entire experiment. 

Following van der Waals* principle of similarity, many thermophysical properties 
for several species seem to be correlated in the sense that one can introduce suit- 
able scale factors in the dependent and independent variables so that properties 
so reduced are congruent to some universal curve (i.e., become universal). This 
applies to normal as well as near critical behavior. But since flow and heat 
transfer behavior depend on these properties, it seems reasonable to expect that 
for similar geometric arrangements it should be possible to describe flow and heat 
transf?.- in such a way as to be, more or less, congruent to come universal curves. 

In general then, one would like to develop solutions which are Independent of BOTH 
the physical nature of the coordinates and the thermodynamic character of the work- 
ing fluid. As an attempt toward accomplishing this, we shall first review some 
similarity principles for the thermophysical properties of fluids in section A, 
with emphasis on the critical region in section B and cite some related computa- 
tional programs in section C,* second, we shall provide some examples of physical 
processes which have been related using similarity principles, in section D; and 
third, in section E we use these similarity principles to demonstrate that two- 
phase choked flows for a simple fluid a can be computed using a simple fluid p 
independent of the thermodynamic character of the other fluids a and p thereby 
demonstrating the process independence of BOTH coordinates and working fluid; our 
conclusions arc summarized in section F. 


A. THE PRINCIPLE OF SIMILARITY IN DESCRIBING 
THERMOPHYSICAL PROPERTIES OF FLUIDS 


Al. Background 

Most engineering texts on therraodynamieg treat the subject of corresponding states 
to some degree for both single and multicomponent fluids, while some consider a 
possible formulation of corresponding states for transport properties. However, 
the actual practice of predicting such fluid properties as a function of the 
thermodynamic state variables requires a great deal of computing capability and 
accurate PVTXqj^ data for developing fundamental equations. In this section, we 
will mention some approaches. 
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A2. Tharmophyslcal Propcctles of One~CompQnent: Flutda 
A2.1 Themodynamtc ProporCics 

Ag an example of the traditional corresponding states methods, consider van der 
Waala' equation of state describing fluid a. 


p . -RT ^ JL 
v-b ^2 


In reduced form it becomes. 


R 


R 






(A.l) 


(A. 2) 


vhere 


Z 


c 







(A. 3) 


and in this form it not only describes fluid a but fluid (3 as well. This ap- 
proach is correct for any system of molecules whose potential can be specified by 
two system-dependent parameters only, 


Mollerup (1975 and 1977) proposes to relace the properties of two pure substances 
by the equations : 


G„(V.T) - - RT In h 


(A, 4) 
(A. 5) 


where o denotes the reference fluid and the properties are termed configurational 
because they evolve from the configurational part of the partition function and are 
conformal in that they have similar molecular potentials and obey the principle of 
corresponding states; for example, both are adequately described by a two parameter 
Lennard-Jones type potential. The functions f and h are dependent on the ratio 
of the critical temperature and volmne 


f = 




(A. 6) 


while 9 and (I> are shape factor functions that depend on the eccentric factors 
of fluid a and fluid o and that are weakly dependent on temperature and volume. 
(See Leach, 1966, 1967; and Eqs. (E.5) and (E.6)). 


Mollerup (1978) has used the Ping-Roblnson equation to predict properties of the 
quantum fluids hydrogen and helium using accentric factors of approximately -0.2187 
and -0.3897 for hydrogen and helium respectively along with the extended corre- 


iWote that the related forms become 


P 

o 




ph; T 


T 

f 


5 


sponding states theory, 
by 


Mollerup finds that the shape factor can be represented 

(A. 7) 


1 + m +(m - in) l/^R 

a o a • g 


1 + m 


for quantum fluids, where is a quadratic function in eccentric factor 

He has verified this representation at reduced temperatures up to Tp < 26 for hy- 
drogen and Tg < 40 for helium. 


A2.2 Transport 

Dimensional analysis^ and corresponding states can be combined to determine suit- 
able scale factors for transport coefficients. 


1, Viscosity data for many fluids can be grouped over a large range in experiment- 
al conditions (Dean and Stiel, 1965; Jossl Stiel, and Thodos, 1962; Stiel and 
Thodos, 1962) using the parameters 

(n “ ~ (A, 8) 


where 




,2/3j,l/2 


(A. 9) 


as illustrated in Fig. 1. 

2. Thermal conductivity data were grouped using a similar parameter (Misic and 
Thodos, 1961; Schaefer and Thodos, 1959) 


^ “ ^0,0^ ^c^2 


(A. 10) 


where 

c 

3. Binary diffusion coefficients were grouped using the parameter (Stattery and 
Bird, 1938) 


where 




(A. 12) 


pV3„l/2 

^3 “ °^/6 

c 


(A. 13) 


2Dynamically similar processes based on the observation that physical laws are in- 
dependent of the particular system of units employed. 


6 


4. Surface tension data con be grouped with the parameter (Brock and Bird, 1955) 


^s^4 


£(\) 


where 



p2/3^1/3 


(A. 14) 


(A. 15) 


Hanley (1977) and Rainwater (1977) applied these configuration functions to thermal 
conductivity and viscosity and proposed that the thermal conductivity and viscosity 
of pure components could be related as^ 

where 



The use of configurational functions f and h along with an accurate reference 
equation of state greatly enhances the users capability to accurately relate 
thermophysical properties, it does however have the disadvantage of being more 
complex and one needs to determine the configuration factors. 

More recent representations of the viscosity surface (in terms of density virlals 
with temperature dependent coefficients) are given by Hanley, McCarty and Haynes 
(1974a) using a corresponding states approach to correlate the results for sever- 
al simple fluids. The procedure is more accurate but can not be applied to as 
many different fluids as some more conventional principles which have been success- 
fully applied not only to simple fluids but to long chain hydrocarbons and polar 
molecules (Dean and Stiel, 1965j Jossi, Stlel, and Thodos, 1962; Stiel and Thodos, 
1962) , 


A3. Thermophvsical Properties of Multicomponent Fluids (Fluid Mixtures) 
A3.1 Perfect Gases 


Gordon and McBride (1971), Zeleznik ani Gordon (1968), and Svhela (1962) have 
formulated a program which calculates the chemical equilibrium composition and mix- 
ing properties of multicomponent fluids. The program can handle fuels and oxidants 
constructed from 15 different chemical elements. There are some 600 different pro- 
ducts of combustion whose data are stored in the computer program, and any given 


3using the related forms of Eq. (A.4)-footnote, we can relate traditional parame- 
ters and the extended corresponding states parameters 
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computation con handle any 200 o£ them. The program includes solid and liquid con- 
densed phases ond ionization, Bittker and Scullin (1972) have a code for reaction 
kinetics. These programs represent a major contribution to combustion computations 
and are based on free energy functions. The thermodynamic and transport property 
contributions of Kestln (1979) and collaborators in formulating accurate represen- 
tations of the thermodynamic and transport properties of pure component fluids and 
fluid mixtures are well known. Most recently they have developed a formula for 
the transport properties of natural gas (Abe and Kestin, 1979). 


A3. 2 Extended Principle-Multicomponent Fluids 
A3. 2.1 'fliermodynamic Properties 

Mollerup (1975, 1977, 1978) used an extended corresponding states principle to 
correlate and predict properties of some fluid mixtures with considerable accuracy. 
In this procedure, a reference substance is chosen, usually the most dominant of 
the mixture with the most accurate PVT surface description, 

f) (A. 19) 

*a> ■ S„(k.|) + r/e \ 1" »a ‘ 1" >’) 

The form is similar to that for pure substances except that wc need to account for 
the energy of mixing and introduce some "mixing rules" which for our case are van 
der Waals' two component rules. While there are certainly other potentially appli- 
cable versions, these appear to be the simplest and most effective. We must now 
redefine the functions f and h: 


""EE 

a B 


(A. 21) 


where 


fh =■ V V X x^h „f „ 
/—J a |3 0.8 aB 
a 3 


(A. 22) 




rh^*^^ + h 


^ap “ 


acL 


1/2 



(A. 23) 
(A. 24) 


and where and are binary interaction parameters. Wliile the procedure 
is complex, Mollerup has already determined many interaction parameters of practi- 
cal interest, for LNG, and has written a computer program to p'redict these proper- 
ties . 
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A3. 2, 2 Transport Properties 

Honley and Rainwater (1977) applies these configuration functions and proposed 
that thermal conductivity and viscosity of mixtures could be related in a manner 
similar to pure components, Eqs. (A. 16) and (A. 17), 

A(p,T,x^) = (A. 25) 

t1(P,T,x^) = q^(ph,|)F^ (A. 26) 

where f and h ore defined as above, but Fr| are redefined.^ For the 

thermal conductivity, Hnnloy proposed an expression based on tlie theory of Enakog 





(A. 27) 




°a,0 


(A. 28) 


where 


^a,0 


(bp) 

Tb^ 


jji-l/exp(p^/p)^ 


■'a,0 ) 

o 


2 “ + 1.2 + 0.775 bpx 
k’I s ^ + 0.8 + 0.761 bpX 


F^, now can be written as 




Cl ■) / 


2/3 


a, 0 




(A. 29) 

(A. 30) 

(A. 31) 
(A. 32) 

(A. 33) 
(A. 34) 


Alt should follow that Frondtl number could be related using these principles, 

Pr = mass diffusion/ thermal diffusion = v/a (for const, properties, Pr= Cp(r|/A)) 

Pr(p,T)^ = Pr(ph,T/f)Q 

or 


[Pr(p,T) - Pr(T)]^ = (dPr)^ 



B. SIMILARITY OF 'rHERMOPHYSICAL PROPERTIES OF FLUIDS 
IN THE CRITICAL REGION 
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B . I Background 

Many thermophysical properties exlilblt anomalous behavior in the vicinity of the 
critical point. For instance, the .isothermal compressibility, the thermal expan- 
sion coefficient and the specific heat of fluids all diverge at the critical point, 
critical enhancement effects arc also encountered in the behavior of the thermal 
conductivity and the viscosity near the critical point. 

It has always been difficult to formulate correlating equations that accurately in- 
corporate the singular behavior of the various thermophysical properties near the 
critical point. However, in recent years considerable progress has been made In 
the art of developing representative equations for the thermodynamic and transport 
properties of fluids in the critical region. 

The critical point is a point of morginal thermodynamic stability. In the vicinity 
of the critical point, large-scale density fluctuations arc present in the fluid. 
The spatial extent of these fluctuations may be characterized by a coherence or 
correlation length Near the critical point, this correlation length becomes 

much larger than the Intermolocular interaction range. One can Imagine that the 
critical enhancement effects arc caused by the static and dynamic behavior of a 
system of clusters. Different states in the vicinity of the critical point corre- 
spond to different cluster sizes and the thermodynamic states can be interrelated 
by scaling laws. Moreover, the properties of a system of clusters become indepen- 
dent of the nature of the individual molecules. This principle is currently known 
as critical-point universality (A. Sengers and j, V, Sengers, 1977). 

Recent developments have demonstrated that these principles become rigorously valid 
only in an asymptotic small region around the critical point, a region too small to 
be of interest in engineering applications. However, using some seml-emplrical 
adaptations, the approach can nevertheless be fruitfully used to formulate repre- 
sentations that can be used In practice in a larger range around the critical 
point (J, V. Sengers and J. M. H. Sengerc, 1978). In this section we consider 
these somewhat empirical generalizations of the principles of scaling and univer- 
sality for representing thermophysical properties of fluids in the critical region. 


B2. Thermophysical Properties 
B2 . 1 Power Laws 


Let F be the Helmlaoltz free energy, S the entropy, ^ the chemical potential per 
unit mass and Cy the heat capacity at constant volume. We also find it conve- 
nient to introduce a symmetrized isothermal compressibility defined as 


X 


T 




(B.l) 


The thermodynamic properties are made dimensionless by expressing them in terms of 
the critical temperature I'c, the critical density P(,, and the critical pressure 
Pc’ 
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T* » T/T , p* » p/p^, P* - P/P^, F* - F/P.V, 

C V C C 


“ MP^/P^, X* = Vc^P^ “ ST^P^V, C; » C^T^/P^V 


(B.2) 


Note that the reduced extensive properties F*, S* and C* are all taken per unit 
volume rather than unit mass. The reason Is that the si;^jular part of the exten- 
sive thermodynamic properties per unit volume appear to be symmetric functions 
(Vincent, .nl-Mlssonl, Sengers, and Green, 1969; J.MtH. Songers, 1974) of (p - p^). 
As a consequence, Xt le also a symmetric function of p - Pq, while the compress- 
blllty k-p s p"1-(3p/5p)x V. Sengers, 1973). We also define the re- 

duced dtfferencen 


At* o (t - T )/T 
c c 

= Cp - Pc^/Pc 


> 


£ii* - [P(P»T) - P(Pc.T)]pj,/PgJ 


(B.3) 


To represent the singular thetmodynamlc behavior of fluids in the vicinity of the 
critical point, one first defines critical power laws (Fisher, 1967), Tlie expo- 
nents of these power laws depend on the property considered and the path along 
which the critical point is approached. The special paths in the AT* - Ap* plane 
commonly considered are the coexistence curve or phase boundary (CXC), the critical 
isoebore (CIC) Ap = 0 and tlie critical isotherm (GIT) AT* « 0, Tliese paths are 
Bchematicnlly Indicated in Fig. 2. The critical power laws defined along these 
paths are summarized In Table I, The critical CKponents of the different thermo- 
dynamic properties arc not independent, The scaling laws to be discussed below 
imply the exponent relations 


a « 2 - 0(6 + 1) 

r » 0(6 - 1 ) 


(B.4) 


2.2 Scaling Laws 

In order to fonrulate scaling laws for the critical region, we introduce a scaling 
varioble x defined as 


X = AT*/|Ap*|^'^^ (B.5) 

This scaling variable is chosen such that the coexistence boundary corresponds to 
a constant va3.ue of the scaling variable x: 


^CXC 


■X , 
o’ 


X 

o 


B-l/S 


Instead of x, it is slightly more convenient to use the variable 


X + X 

o 



(B.6) 


CB.7) 


II 


Curves of constant y iire schema ttcnlly indicated in Fig. 3. The coexistence 
curve corresponds to y 0, the critical isotherm to y « I and the criticol iso- 
chorc to y >» ». 

In the critical region the thermodynamic properties arc considered as a sum of two 
types of contributions, namely one typo of contribution then represents the singu- 
lar behavior at the critical point and sotisfJ.es scaling laws and another type of 
contribution that it. a regular or "background" contribution. The sealed expres- 
sions for a number of thermodynamic functions ore summarized in Table 2 (J. V. 
Sengers and J, M, H. Sengers, 1978). Here the slngulor behovior of the thormo- 
dynomic properties is expressed in terms of two scaling functions h(y) and a(y) 
which are related by the differential equotion (J. M. H. Songeri, Greer, and J. V. 
Sengers, 1976) 


?h(y) = + l)o(y) + (1 - y) (B.8) 

The function h(y) is normalized such that h(l) = 1. 

As an example, we consider the scaling low for the chemical potential 


V' “ V'[A:i*|^"^mi(y) (B.9) 

Normally, the chemical potential Is a function of two variables, density end tem- 
perature. The scaling law (B.9) Implies that when the data for Ap* are divided 
by the data should collapse onto a single curve as a function of the 

scaling variable y. In Fig. A we show such a scaled plot for steam (J. M. H. 
Sengers, Greer, and J. V, Sengers, 1976). 


To use the scaling lavi?8 in practice, one needs an explicit expression for the seal 
ing function h(y). An expression proposed for this purpose is the so-called NBS- 
or MLSG-equation (Sengers, Grear, and, Sengers, 1976; Sengers and Sengers, 1977; 
Vinccntinl-Mlssoni, Sengers, and Green, 1969) 


h(y) 


\ 1 + E 




(r-i)/23 


(B.IO) 


where E is a constant to be discussed below, Tlie advantage of this expression 
is that it enables one to calculate the compressibility directly as a func- 

tion of density and temperature in the critical region. The disadvantage is that 
it cannot be analytically integrated to yield an explicit expression for a(y) and 
hence, for the pressure. This problem is solved by introducing two auxiliary 
("parametric") variables, replacing the variables AT* and Ap*. For a discus- 
sion of the parametric equations for the thermodynamic properties in the critical 
region, the reader is referred to a review published elsewhere (Sengers and 
Sengers, 1978). 


It follows from (B.5) and (B.9), that all critical exponents for the power laws in 
Table 1 are related to the exponents 3 uud 6. Tliese relations were earlier 
given in (B.4). Tlie amplitudes B, P and P' of the critical power lows for the 
coexistence curve and the compressibility arc related to D and x^ by 


o 


(B.lla) 
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r « 


X / 

f(- 


Li_E^(r-i)/20 

E J 


(B.Ub; 


r' 




(1 + E) 


(r*i)/2|3 


(U.llc) 


2.3 Unlveraolttv 


The principle of universality asserts that the critical exponents end the scaling 
function h(y) should be universal, that is, the same for all 3 -dimensional systems 
with short range forces. Tlic universal values for the critical exponents of this 
universality class have been calculated theoretically (Baker, Nickel, and Mciron, 
1978; Guillon ond Zinn-Justin, 1977). The only system dependent constants in the 
scaling laws arc the coefficients D and Xq which specify the amplitude of the 
critical isotherm and the coexistence curve, respectively. 

It turns out that the principle of universality in its strict sense applies in s 
small region AT* < 10"^ near the critical point which is insignificontly small for 
engineering applications (Sengers and Sengers, 1977). However, an empirical exten- 
sion of the i'f’<versality principle as applied to fluids near the critical point is 
obtained b" *.n‘ ''oducing effective fluid exponents and an effective value for the 
constant iu (B.IO) (Sengers, Greer, and Sengers, 1976; Sengers and Sengers, 
1976; Sfmgers ond Sengers, 1978) 

a « 0.100, 3 «= 0.355, r = 1.19, 6 = 4.352, E » 0.287 (B.12) 

As an illuEtrotlon of this practical application of universality, we show in Fig, 5 
scaled chemical potential data (Sengers, Greer, and Sengers, 1976; Sengers and 
Sengers, 1978) for %e, ^He, Xe, CO 2 and H 2 O. The data can opproxiraately be re- 
presented by 0 single function of the scaling variable y = (x + Xq)/xq. The data 
cover a range in temperature and density corresponding to 

SXIO"'^ < jAT*| < 3X10"^, |Ap*| < 0.25 (B.13) 

The principle of corresponding states would imply identical behovior when the 
thermodynamic properties are reduced using the critical parameters Tg, Pc, Pg. 

The fluids Included in Fig. 5 dn not satisfy corresponding states. However, the 
universality hypothesis is a generalizotlon of corresponding states, saying that 
universal critical behavior is obtained using only two additional constants, U and 
Xo, besides the critical parameters themselves. In Table 3, taken from an ear- 
lier review (Sengers and Sengers, 1978), we have listed these parameters for a 
large number of fluids . 

The principle of universality applies to the compressibility Xm and to the singu- 
lar contributions of the other thermodynamic properties (see Table 2). It may be 
noted that (5P/dT)p is finite at the critical point, while the compressibility 
Xx diverges strongly. It has been found In some practical applications that al- 
ready a significant improvement is obtained if Xx it! calculated from the scaling 
laws even if the equation of state and, hence (Sp/3T)p are represented by more 
conventlonul equations. 
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3. Correlation ^ ncnath 

Ai\ important quantity in die critical region is the correlation length | which 
moaflurea the size of the cluatere. In practice this correlation length may be 
estimated from (Sengers and Sengnrs, 1978) 

I = (B.14) 


where T is given by (D.llb). The critical exponent v for the correlation 
length in related to the thermodynamic exponents p and y* by the hyperscaling 
relation 


V a (23 r)/3 


(B.15) 


Using the values adopted in (B.12) one finds v » 0.633. The scale factor So for 
the correlation length can be related to the critical power amplitudes B and T 
by (Sengers and Sengers, 1978) 



(B.16) 


where is a universal constant whose value is approximately equal to 0.7. 

Values of the corielation length scale factor Iq for a number of fluids are also 
Included in Table 3. 


4. Transport Properties 


4.1 Thermal Conductivity 

The thermal conductivity of fluids shows a pronounced critical enhancement In a 
large range of temperatures and densities around the critical point (Sengers, 1972, 
1973). To .iCC''";rt for this behavior the total thermal conductivity is decom- 
posed as the sum of a normal part A and a critical part 

A = A + M (B . 17 ) 

The idea is that the normal part A = A + ^0,0 represents the thermal conductivity 
without the extra contributions due to the critical fluctuations. This part, 
therefore, can be represented by customary correlating procedures, section sec- 
tion A. 


The theory of dynamical critical phenomena predicts that the critical enhancement 
M will diverge asymptotically at the critical point as (Hohenberg and Halperin, 
1967; Kawasaki, 1970; Sengers, 1973) 


M - AA 

c 


67nil p2 VaT/p ^ " 67m& p*2 \ST*/p ^T 


(B.18) 


where kg is Boltzmann's constant, q the shear viscosity and A a universal 
constant of order unity. Values predicted theoretically for this constant A span 
the range 1.0 < A < 1.2. 

An empirical generalization of (B.18) to cover the entire range of temperatures and 
densities of the thermal conductivity enhancement, is obtained by writing (Hanley, 
McCarty, and Sengers, 1974; Hanley, Sengers, and Ely, 1976), 
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» A\^fj^[AT*,Ap*| (B.19) 

where is given by (B.18) and where f(AT*,^*) la an empirical damping func- 

tion auch that 


lim f. (AT*,Ap*) « 1 
AT*-*0,Ap*-.0 ^ 


and 

lim f (AT*,Ap*) = 0 
AT*-Ko,Ap*-«o ^ 

An exf^mple of a suitable damping function is (Hanley, McCarty, and Sengers, 1974; 
Hanley, Sengers, and Ely, 1976) 

f^(AT*,Ap*) = exp [c^(AT*)^ + Cp(.!ip*)^]J (B.20) 

with = 18.66 and Cp = 4.25. Equations of the type given in (B.19) have been 
used to represent the thermal conductivity of a variety of fluids (Hanley and 
colleagues, 1974a, 1974b, 1976, 1977a, 1977b) including steam (Basu and Sengers, 
1979). 


4.2 Viscosity 


Theory predicts that the shear viscosity q of fluids also diverges at the criti- 
cal point (Hohenberg and Halperin, 1967; Ohta and Kawasaki, 1976; Perl and Ferrell, 
1972). Hence, we write in analogy to (B.17) 

q = q + Aq (B.21) 

where q = q + Hq 0 normal viscosity and Aq the critical viscosity en- 

hancement. It turns out that the critical viscosity enhancement Is a multiplica- 
tive enliancement, that is, Aq is itself proportional to q (Ohta, 1977). ^There- 
fore, in practice, we prefer to consider the viscosity ratio Aq/q or q/q (Basu 
and Sengers, 1979). 

Unlike the critical thermal conductivity enhancement, the critical viscosity en- 
hancement is weak and limited to a small range of temperatures and densities around 
the critical point (Basu and Sengers, 1979). Hence, the enhancement may be ne- 
glected for many applications as »'as done in the viscosity equation recently 
adopted by ZAPS for steam (Hendricks, McClintock, and Silvestri, 1977; Nagashima, 
1977; Schleffer and colleagues, 1978). 

A more detailed discussion of the critical viscosity enhancement of fluids near the 
gas-liquid critical point is presented elsewhere (Basu and Sengers, 1979; Basu, 
Sengers, and Watson, unpublished). It can be estimated as 

qi§ > i"l 

> (B.22) 

qil < M 

where cp^ is a universal critical exponent and a parameter depending on the 
fluid. The effect of the critical enhancement is small and accurate knowledge of 


q = q(q^l)'^’^ for 

q = q for 
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Cpj^ and qj^ is not needed for practical applications. As a guide we recommended 
an effective exponent value (Basu and Sengers, 1979) 

cp^ = 8/157T^ e 0.054 (B.23) 

while the parameter qi may be estimated from (Basu and Sengers, 1979; Oxtoby and 
Gelbart, 1974) 




V 

SttbAp 


^ Vc _r /Sp*V 

Ut; ^2 " ,2 VW 




>: 


while, as "a rule of thumb," (3p*/3T*)p 6, see Eq. (D.4). 


(B.24) 


C. NUMERICAL PLOTTING CODES 

Before preceding, we must caution the rejder the’, accurate representations of the 
PVT surfaces ore a requirement, for a meaningfil transformation of any aalution 
to physical reality. 

Equations of state (EOS) are available from NBS, lUPAC, HTFS, NASA, Tech. Unlv. 
Munich, etc. and mixtures from Mollerup Tech. Unlv, of Denmark, AGA, NBS, NASA and 
the Petro Chem. Ind., API.^ 

With the advent of the computer and plotting programs, such as those developed by 
Trivisonno and Hendricks (1977), and the thermophysical properties routines for 
several fluids (Hendricks and colleagues, 1973, 1975), or any of those supplied 
from the organizations list above, thermophysical properties may be accurately and 
quickly compared on such complex charts as T-S, H-S, H-P, F-T, G-T, etc., see 
Fig. 6. These charts, which are useful in cycle analysis, also serve to rapidly 
assess descriptions of the PVT surface and its derivatives which, in turn, can 
establish corresponding states reduction parameters as required. The recent work 
of Mollerup extends the input parameter options of his mixtures program to (P, cp, 
Xq^) where tp = T, V, H, S and (H, S, Xq^) . 


D. EXAMPLES OF THERMOPHYSICAL PROPERTIES RELATED 
TO DYNAMIC PROCESSES 

Dl. Thermophysical Properties 

Dl.l A classic example of a thermodynaidc property related to a dynamic process is 
sonic velocity. The dynamic wave motion is described completely and entirely by 
curvature of the PVT surface representing the fluid through which the wave is 
passing. 


c 


2 


s 



(D.l) 


Thus sonic velocity is fixed by state and not process. This concept is modified 
somewhat at high frequencies. 


5nBS ~ Gaithersburg, MD and Boulder, CO; lUPAC - Imperial College, I,ondon; HTFS - 
AERE Harwell, England; NASA - COMSTC U. of Georgia, Athens, GA; AGA - Arlington, 
VA; API - Washington, D.C. 
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D1.2 JouIc-'Thomson coefficients arc also thermodynamic properties which are related 
to on expansion process. For an Iscnthalplc expansion, 


D1.3 Hendricks, Baron, and peLler (1972) demonstrated that for engineering purposes 
the reduced Joule-Tliomson coefficient, see Fig, 7, was useful in grouping results 
for various simple fluids 





OR 



where 




o,c 



(D.3) 


(D,4) 


Figure 7 illustrates a reduced Joule-Thompson coefficient plot, 

D1.4 Leinhard and Korimi (1977) and Alamgir applied the corresponding states prin- 
ciple to I an der Waals spinodal locus (Leinhard and Karirai, 1977) and surface ten- 
sion (Leinhard and Alamgir) for several fluids, see Fig, 8, 

Dl.5 Simoneau and Hendricks (1977) used the principle to extend the calculations 
for equilibrium fluid properties into the raetas table regions to the equivalent van 
der Waals minimum pressure locus or vapor spinodal, for use in two phase choked 
flow computations, see Fig, 9, 

Dl,6 Parameters such os (h In p/Sepj) = 3j where cpj = T, H, E, A, or can be re- 
duced to a single plot of sufficient accuracy for most heat transfer calculations. 
This has been demonstrated by Hendricks (1979). Figure 10 represents a reduced 
volumetric expansion factor plot = 3l/(ZcTc)"^. These few examples serve to 
illustrate some of the principles of similarity. 


D2. Heat and Mass Transfer Processes Related throuRh the Principle 
D.2,1 Choked Flows 


Hendricks (1973, 1974) has shown that two-phase choked flow rates of several simple 
fluids can be grouped using the reduced mass flux (see Fig. 11(a)) 


where 


W. 


(D.5) 


6lt should follow that from Eq. (A.4)-footnote the extended corresponding states 
relation is: 

a (p,T) = o (ph,T/f)(f/h^''^) 


17 


W* 



( 1 CLosslcal Fluids 
f(Tj^) Quantum Fluids 


(D.6) 


also in good agreement are the reduced pressure profiles of Fig. 11(b) where the 
ordinate represents the ratio of the nozzle throat pressure to that of an isentro- 
pic gas expansion. 

The parameter has been useful in grouping and predicting moss flux for sever- 
al geometries in addition to multiple fluids, see Simoneau and Hendricks (1979). 

Hendricks (1979a, 1979b) and Poolos (1979) have used this normalization procedure 
to group data for free jet flows through tubes with Borda and sharp edge inlets. 
For (Hendricks and Poolos, 1979) the tube was short, L/D « 2, but four inlet geom- 
etries were used; circular, square, rectangular and triangular. 

More recently, Hendricks (1979c) has shown that for gases flowing through a simu- 
lated shaft seal configuration, a modified form 



(D.7) 


reduces to a universal curve for simple gases; for this simulated shaft seal geome 
try, c is nearly constant at 0.16. For a nozzle, c = 1/5. 

In convergent tapered passages, the mass flux rate is reduced using the area 

^1 ' ^^^exit ^inlet^'^^ 

W = rii/Aj^ (D.9) 

Mich this modification, the above relations are valid for convergent tapered pas- 
sages as well. 


D2.2 Heat Transfer 


Hendricks (1977) has shown that there is good evideni.e that the transposed critical 
temperature, used in near critical heat transfer problems, can be related 
through the corresponding states principles, 


7a more general form is ; 

w* . wJv/fTh 

and for the ideal gase case, 

= Const. I f ^ ) 

a R \k + 1/ 

where k is the isentropic expansion coefficient. 

®T* represents the loci! of maximum specific heats 


(k+l/k-1) 
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The transposed crltlc^il temperature locus has been used by several Investiga- 
tors to classify heat transfer data where behavior was more gas like T > T* or 
more liquid like T < and processes which crossed the locus could possess 

either a maximum or a minimum heat transport. This phenomenon Is more fully dis- 
cussed by Hendricks and colleagues (1962, 1966, 1970). The point Is, that T* ap- 
pears to be a corresponding states parameter, some heat transfer data can be clas- 
sified using this criteria and consequently the grouping of heat transfer and hy- 
drodynamic data in this manner could also follow a similarity principle. 

Two component models are commonly used to represent heat transfer and hydrodynamic 
data where distinct phase transitions occur (e.g., liquid-vapor). However, two 
component models are also used to represent Hel and Hell transitions below the 
temperature (Lee and colleagues, 1978), and In grouping data of classical fluids 
in the subcrltlcal and supercritical pressure states (Hendricks and colleagues, 
1962, 1970), Similarities in wall temperature profiles and flow stability between 
the subcrltlcal and supercritical pressure heat transfer data for forced convec- 
tion and the wall temperature-heat flux profiles and stability in "pool boiling" 
have been noted by several investigators (Graham, Hendricks, and Ehlera, 1964 j 
Hendricks and colleagues, 1962; Sabersky and Hauptmann, 1967). 


D2.2.2 Subcrltlcal Convection 


Hendricks and Papell (1977) used the principle to group forced convection sub- 
cooled nucleate boiling data for neon, nitrogen and water. The reduced tempera- 
ture difference was expressed as 


where 


and 


AT„ = 


- f) 


1/4 


4p, 


R 3 


qr, exp 


R 






OJ = -1 - log P, 


V 


T^.0.7 


< ■ '’tsCs" 4 >=tg 


:) 


1/4 


is the Kutateladze expression for critical heat flux 


(D.IO) 


(D..11) 

(D.12) 


(D.13) 


D2.2.3 Volumetric Expansion 

Yaskiii and colleagues (1977) used the volumetric expansion parameter In the form 
of (1 + At) in. successfully grouping heat transfer data for fluid helium in- 
cluding the near critical region. A similar parameter [l + | (3 In p)/(SH)|p Ah] 
was used by Hendricks (1977), to group fluid nitrogen and fluid oxygen data. How- 
ever it did not group the data as well as the work of Yaskin and colleagues (1977). 
Here AT ^ - Tg and AH = - Hu. 
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Hendricks (1977) showed that heat transfer data for fluids nitrogen and oxygen 
wore qualltotivoly related using the principle; however, assessment of the effects 
of buoyancy, relative velocity porameters, etc, which also grouped these data, 
must be mode before the simulation can be made quantitative. 

It was shown by Simoneau and Hendricks (1964), that for gases, the heat transfer 
coefficient can be written 


hj. « CG°'® d’°'^ 

and it con be shown that for an ideal gas 

[l + Bi<T„ - - vV^ 


(Note: C related to A(*) This work was reviewed (Hendricks, 197 9d), and many 
sets of heat transfer data can be grouped using volumetric expansion parameters 
(pjj and the reduced Nusselt number Nuj^, see Fig. 12: 


Nu 

^ R Nu , ^k 
calc 



(D.16) 


where 


Nu , = 0,023 Re®’® pr®’''^ (D.17) 

calc 

= 1 + = 1 + 3|^(tp^^ - tp^)^ no sum on k 


cpj^ = T fp'2 = F cp^ = E ■) 

tp2 = H cp^ = G 'Pg = J 


(D.18) 

(D.19) 


Using the Gibbs energy parameter (p^, the data group well indicating that if the 
process were constructed using Gibbs potential, a better understanding would ensue. 
Recall that the scaled equation of state in the critical region is based on chemi- 
cal potential which reduces to Gibbs potential for pure fluids. 


Tlie coefficient, p^, has been shown (Dl.6), to obey the principle of corresponding 
states; for k = 1 the normalizing parameter (Zj.Tu)"^, as noted before, is suit- 
able for engineering calculations. 


D2.2.4 Volumetric Expansion and Density Fluctuations 

Hendricks and colleagues (1979) have demonstrated that the parameter pj is re- 
lated to density fluctuations and in turn to turbulent flows with substantial tem- 
perature and or pressure gradients. Further, the fluctuation components are re- 
lated through the equations of state where the density can be written 
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- ^ = O^h' + k^.P’ + (D*21) 

where the "primed" terms represent the fluctuating ten.’»'J. JCntroduclng fluctuation 
density (p => p + p'), pressure (P » P + P'), enthalpy (}• ® H + H*), and velocity 
(u ■= u + u') Into the conservation equations cited In the INTRODUCTION, the turbu- 
lent momentum and energy transport, for developed pipe flow, were written 
(Hendricks and colleagues, 1979), 


-Tj. “ p u'v' + V p'u' + p'u'v' 

(D.22) 

"«!(. “ P v'H' + V p'H' + p'v'H' 

(D.23) 


At very small temperature differences, die convective heat transfer is enhanced and 
It follows that the Rayleigh number would be an Important parameter 

B pCppiVT 

See for example (Hendricks, Simoneau, and Smith, 1970) and the careful work of 
Irie and colleagues (1978), which would follow more closely the scaling laws; com- 
pare also the work of Lee and colleagues (1978) and Kllpplng, 


E. POSSIBLE GENERALIZATIONS 
El. Background 

The usefulness of the similarity principles in grouping results for processes with 
heat and momentum transport has been demonstrated in the previous sections. The 
fact that the results obtained with cryogenic fluids can also be applied to a wide- 
ly different fluid such as water is extremely useful. 

Following the extended corresponding states principles of sections A and B, in this 
section we propose to "scale" fluids as 

% ■ I”,’' ■'o ■ \ " ''o “ 

and transform thermophysioal properties as 

cp"^ = /y(h,f,®,0.F^,F^— )cp 

such that solutions of the fundamental equations of motion could be transformed in- 
dependent of fluid characteristics.^ 


9To assign a characteristic length common to both molecular and continuum mechan- 
ics, it will be necessary to separate the effects of h, F^, F^ (Eqs. (A.4)- 
(A. 16-17). Then perhaps one could assign 

L = L (e /CM 
a o\a/^o/ g 

however the assessment of is difficult as transport property measurements 

themselves represent solutions of the conservation equations and the character 
of fg must be such as to give Lq_ independent of cluster size where appropri- 
ate; its definition is beyond the scope of this paper. 
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E2. Calculation of Critical Moss Flux Ualng the Goneralized Transform 

The thermophysical properties program o£ Hendricks and colleagues (1975) • which 
was based on the work o£ Binder (1970), was used to both compute the properties o£ 
the re£erence £luid and tho necessory constants £or the transiorm. The anolytical 
program used by Simoneau and Hendricks (1979) was adapted to calculate the flow 
rates ond pressure ratios for two phase choked flow through a nozzle. The pres- 
sure and temperature of the working fluid were converted to the pressure and tem- 
perature of the reference fluid; computations then proceeded as usual for a single 
co.jponent reference fluid. 


The conversion of pressure and temperature follows 


where 


a 


R 


r 

CL f 

P B P — = 

a o „c h 
*o 


h 


‘R 


a 


T « fX 
a o 


T 

£ = — 0 


0 = 1 + (CO^ - 0J^)[b^(Vj^ - b^) - b 3 (Vj^ - b^) In T^] 
e = 1 + (to^ - Cb^)[a^ - a^ In Tj^ + (s 3 - a^/\) (Vj^ - a^)] 


(E.3) 


(E.4) 


(E.5) 

(E.6) 

(E.7) 


The values of <J> and 0 were adopted from McCarty (1977). 

Some calculated reduced critical mass flux values are given in Table 4 for oxygen 
or nitrogen flowing through o nozzle. The computed values of reduced mass flux are 
in good agreement, as can be seen in Tables 5 and 6 which were constructed using 
the various reference fluids methane, nitrogen, oxygen. Tlie calculated values are 
compared to the experimental data. Fig. 12, of Hendricks, Simoneau, and Barrows 
(1976) in Figs. 13 to 19. In general the nitrogen and oxygen data points are in 
excellent agreement, independent of the working fluid, providing they are normal- 
ized (Hendricks, Simoneau, and Barrows, 1976), see Fig. 12. The theoretical calcu- 
lations are for the most part, in good agreement independent of the working fluid 
and the reference fluid used to calculate the results. Tlie agreement between oxy- 
gen d/ita and theory using oxygen properties is in general quite good for this com- 
plex flow and the same comment can be made for nitrogen. 


These results lend credence to the similarity principles delineated earlier and in 
part substantiate their validity. Tlie work for the thermophystcal properties peo- 
ple is extensive as application of the similarity principles requires accurate data 
bases for classes of fluids. In fluid mechanics and heat transfer accurate defini- 
tive experiments will be required. 


lOxhe homogeneous equilibrium fluid model of the nozzle with no slip between phases 
was chosen because of its relative simplicity and agreement with experimental 
results. 
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F. SUMMARY 

The principle of aimtlarity as applied Cu thermophyaical properties with emphasis 
on the critical region and fluid mechanics, has been reviewed and shown, qualita- 
tively, to group the experimental results of several physical processes. Using 
these principles, it is proposed that solutions can be developed which arc indepen- 
dent of BOTH the fluid charaoterlstlca and the physical nature of the coordinates. 
To illustrate this proposal, two phase choked flow of fluid a was computed using 
fluid 0 where a and p were either methane, nitrogen or oxygen; the results 
ate in good agreement with extensive oxygen and nitrogen data. 
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area 

constants 

scaling parameter, Eq, (B.6) 
specific heat, constant pressure 
specific heat, constant volume 
constant 
speed of sound 

coefficients of temperature and density difference terms, 
Eq, (B.20) 

diffusion coefficient 

constants, Eqs. (B.9), (B.IO), and (B.12) 
diameter 

free energy, Helmholtz 

scale function for thermal conductivity 

scale function for viscosity 

scale function, fluid a to reference fluid 

friction factor 

bounding function, Eq. (B.19) 

configurational Gibbs free energy; mixture 

configurational Gibbs free energy, reference fluid 

configurational Gibbs free energy, fluid a 

Gibbs free energy 

scale parameter, Eqs, (A. 30), (A. 27), and (A. 28) 
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H 
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h(y) 

k 
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Pr 
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V 

W = pu 

“* " V^TTh 

X 


body force 
static enthalpy 
totol enthoLpy 

scale function, fluid a to reference fluid 

heat transfer coefficient 

scaling function, Eq. (B.8) 

isentropic exponsion coefficient 

Boltzman's constant 

Isothermal compressibility 

molecular weight 

mass flow rate 

Prandtl number 

pressure 

scale parameter, Eq. (A, 29) 
heat flux 

reference heat flux, Kutateladze, Eq. (D,13) 

fluid scale parameter, Eq. (B.22) 

gas constant 

entropy 

temperature 

transposed critical temperature (pseudo critical tem- 
perature) 

time 

reference time 
velocity components 
specific volume 
mass flux 

flow normalizing parameter 
coordinate scale 


X 


scaling variable, Eq, (B.5) 
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X mole fraction 

a 

y scoling poromoter, Eq, (B.7) 

Z ■ PV/RT comprosalbility 
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0 

K 
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A 
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V 

a 

i 

h 

h 

&3 

^4 

P 

a 

a 

T 
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omplitude of power law relation for specific heat 

Bcaling low parnmetorfl, Eqs. (B.4) and B.12) 

volumetric or energy expansion coefficient 

scoling parnmetei, Eq. (D.llb) 

scaling parameter, Eq, (B.llc) 

viscosity 

shape factor 

Enskog parameter, Eqs, (A. 33) and (A. 34) 

numericol level constant, Eq. (B.19) 

thermal conductivity 

Hel - Hell transition locus 

chemical potential 

Joule-Thomson coefficient 

critical exponent, Eq, (B.14) 

dynamic viscosity 

correlation length, Eq. (B.14) 

viscosity normalization parameter, Eq, (A. 5) 

thermal conductivity .tormalization parameter, Eq. (A. 7) 

diffusion normalization parameter, Eq. (A. 9) 

surface tension normalization parameter, Eq. (A. 11) 

binary interaction parameter 

density 

surface tension 
shear force 
shape factor 
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property pdremeter^ Bqs. (D.19) end (E.2) 
blnory Interaction parameter 
constant, Eq. (E.23) 

scale parameter, Bqs. (A. 27), (A. 28), and (A. 33) 
symmetrized Isothermal compressibility 
correction parameter, Eq, (D.6) 

volumetric or energy expansion parameter, Eq. (D.17) 
accentric factor, Eq. (D.12) 
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RTOUCEO FKISSURI. Rp 

F 19 . I< Compiriioni of choktd How rain at a function of raducad praitura 
tor raducad lamparatura of -<l 7M. 



Fig. IS. rompamons of chokad flow ralas as a function of raducad prtssura 
tor raducad lamparatura of '0. 942. 
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Fig. 16. CoffliMrisons ol chok«) Mow rdts <t i lunclion ol roducod pmturt 
lor rtducod t•mptrltur* ol -1.026. 
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f jg. 17. Comparison ol pressure ratio as a lunclion ol reduced pressure lor 
reduced temperature ol -0. 750. 




RATIO or THROAT PRtSSURi TO STAGNATION 
PRtSSUHt. P,;Pq 



Ttq. 18. Compirison ol prfsturt rtlio ii i function of rtductd prmurt lor 
riducid lomporiturt of ~0. 942. 
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Tig. 19. Compirison of pmsurt ratio as a function ol raducad prassure lor 
raducad tamparatura of -1.02S. 
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